Current-current Fermi-liquid corrections to the superconducting fluctuations on 

conductivity and diamagnetism 
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We analyze the behavior of the superconducting-fluctuations contribution to diamagnetism and 
conductivity in a model system having current-current interactions. We show that in proximity 
to a Mott-insulating phase one recovers an overall suppression of the fluctuating contribution to 
the conductivity with respect to diamagnetism, in close analogy with recent experiments on the 
underdoped phase of cuprate superconductors. 
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It is generally believed that, due to their low superfluid 
densities and short correlation lengths, superconduct- 
ing fluctuations (SCF) in underdoped cuprates should 
be relevant for transport and thermodynamic proper- 
ties. Such SCF have been widely highlighted in the 
pseudogap region by several experimental measurements, 
ranging from diamagnetism^"'* and Nernst effect^'^ to 
paraconductivity^"® . In particular, the survival of a large 
Nernst signal up to temperatures much larger than the 
superconducting transition temperature T^ in the under- 
doped region has been interpreted as the evidence for 
vortex-like phase fluctuations with a Kosterlitz-Thouless 
(KT) character due to the quasi-two-dimensional (2D) 
nature of the system^''^. At the same time, several 
authors^' claimed that paraconductivity in underdoped 
cuprates simply follows the T dependence expected for 
the quasi-2D Aslamazov-Larkin (AL) regime of gaussian 
Ginzburg-Landau (GL) SCF close to Tc^ . This outcome 
motivated various investigations in the attempt to ex- 
plain the experimental data on Nernst signal and dia- 
magnetism within a GL-like framework*''"^^ . 

Regardless of the KT or GL character of fluctuations, 
there are two outcomes of the experiments that are not 
expected for conventional superconductors: (i) the range 
of temperatures where the fluctuation conductivity is ob- 
served does not always match the one where a sizeable 
Nernst signal has been reported^'*, and (ii) the SCF con- 
tribution to the conductivity is about two orders of mag- 
nitude smaller than the fluctuating diamagnetism in the 
same system, as recently pointed out by Bilbro et al. in 
Ref. 13. To be more precise, we recall that in 2D the 
contribution of SCF to the conductivity Sa and diamag- 
netism 5xd can be expressed both within the GL^ and 
KT*^ theory in terms of the superconducting correlation 
length ^(T) as: 



Sxd 






em 



Sa = 



em 



16hd 



so 



(1) 



Here $o is the flux quantum, ^o the low temperature 
correlation length and d is the thickness of the effective 
2D system (i.e. the interlayer spacing for layered systems 



as cuprates). Both ba and bxd diverge as T approaches 
Tc due to the increase of the correlation length ^(T), 
with a T dependence that is power-law within the GL 
approach and exponential within the KT theory. From 
Eq.s (1) one could express the ratio 5a/5xd as 
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where ^^j Cxd ^^^ ^^^ correlations length extracted from 
paraconductivity and diamagnetism measurements, re- 
spectively. Since one would expect that the same lenght 
scale is involved in both cases, Cct/Cv^ = 1: rigth above 
Tc, Sa/6xd should be of order - lO^{nA/T)-\ while ex- 
perimentally it turns out to be two orders of magnitude 
smaller than this*'^^ (see also discussion below Eq. (36)). 
Let us notice that the above discussion holds regardeless 
the nature of the SCF so that the the ratio between 6a 
and 6xd depends only on the properties of the system 
away from Tc- 

In this paper we show that the quantitative disagree- 
ment between the SCF contribution to diamagnetism and 
conductivity could be understood as a consequence of 
current-current interactions in a doped Mott insulator. 
The possible relevance of this kind of interactions to the 
physics of cuprates has been suggested within several con- 
texts, ranging from the gauge-theory formulation for the 
t — J model*'"' to the theoretical approaches emphasiz- 
ing the role of microscopic currents*^'*^. As a paradig- 
matic example we focus on the t — J model within the 
slave-boson approach, where the Hartree-Fock (HF) cor- 
rection of the quasiparticle dispersion leads to a differ- 
ence between the quasiparticle current and velocity in 
the usual Landau Fermi-liquid (FL) language*^. The 
effect of this dichotomy on the GL functional for SCF 
can be accounted for within a general fleld-theory for 
the t-J model*^'*^, which includes fluctuations both in 
the particle-particle (p-p) and particle-hole (p-h) chan- 
nel. By computing the AL contribution to diamagnetism 
and conductivity we find that the role of Landau FL 
corrections differs in the static or dynamic limit. As 
a consequence the SCF contribution to diamagnetism 



and conductivity scales with a different prefactor, lead- 
ing to the suppression of paraconductivity with respect 
to fluctuation diamagnetism when the Mott insulator is 
approached, in analogy with experiments in cuprates. 

Let us start from the slave-boson version of the t ~ J 
model. 



H=-t5Y, iAaCja + H.c.) + jY, S.S, (3) 
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where t is the electron hopping, J is the exchange inter- 
action between electron spins, cl^{cia) is the fermionic 
creation (annihilation) operator and Si = ^J^^i is the 
spin operator, with ^^ = il'^)- The sum is extended 
over all the («, j) nearest-neighbor pairs on a square lat- 
tice, and we use units such that the lattice spacing and 
h = c = e = 1. In Eq. (3) the decomposition of 
the electron operator in a fermionic spinon and bosonic 
holon part has been already carried out^'', and only the 
fermionic degrees of freedom have been retained. Here we 
neglect the boson fluctuations, and we assume that slave 
bosons are always condensed, leading to the suppression 
factor S = 2x/(l + x) of the hopping, scaling with the 
doping X. The interaction term of the Eq.(3) contains 
contributions both from the p-p and the p-h channel. 
We introduce the operators 






^COsk„ 4^q/2^^Ck_q/2,a, (4) 

^sink„c^q/2^^Ck_q/2,., (5) 

ko- 

^7d(k) C_k+q/2,iCk+q/2,t (6) 



where a = x,y and 7d(k) — cosk^, — coskj,, so the inter- 
action term reads 



Y^ss, = -gj: f jE'^^.*^q+*^.'^"0+'^^*'^ 
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(7) 
with g = 3 J/4:. The RHS of Eq. (7) represents the inter- 
action in the p-h density ($'^), p-h current (<&") and p-p 
((f>^) channel, respectively^". 

We decouple (7) by means of the Hubbard- 
Stratonovich (HS) trasformation both in the p-p and in 
the p-h channel. After the integration of the fermions 
the action reads: 



s^j2ij2\^£^\^'\ , |A,|2 
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(8) 
Here cj)'^, (/>* and A are the HS fields, the trace acts 
over momenta, frequencies and spins, k = (k, i£„), q = 
(q, iujm), and £„, ujm are the Matsubara fermion and bo- 
son frequencies, respectively. The Ak.k' matrix is defined 



(in the usual Nambu notation) as 

^k'k = -[ia;„ + sin(is±!s:)„0^ fc_fc,]fo 



e£ - cos(i^)„C k-k' 



where ^k is the bare dispersion 
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The Eq. (9) can be decomposed as A^.k' = ^Gq ^S, 



"fe,fe' 



^k — k' 



Here Gq Sk^k' contains the q = Q saddle-point val- 
/)g, (^g, Aq) of the HS fields obtained by minimization 



of the mean-field action 
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nil: 
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while Sfe_fe' contains the fluctuating parts of the HS 
fields. The standard GL functional above Tc will then 
be given by the expansion of Eq. (8) around the mean- 
field action (11) as Sgl = E„ Tr[GoS]"/n, with Aq = 0. 
As far as the saddle-point values in the p-h channels are 
concerned, one can easily see that </>§== 0, while (/)q 7^ 
satisfies the following self-consistent equation: 
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where ,f{x) is the Fermi function, j3 = 1/T and ^k is the 
quasiparticle dispersion, given by the 73 term of Eq. (9): 

Ck = — (2M + (/'o)(coska; -I- cosky) — /i (13) 

As one can see, the 0g value corresponds thus to the 
standard HE correction of the quasiparticle dispersion. 

In order to compute the contribution Sxil) of SCF 
to the electromagnetic response function we need to in- 
troduce in the effective action also the electromagnetic 
potential A. For the model (3) this can be done via the 
Pierls substitution clci+a -^ clci+ae~^^' ■ which modi- 
fies the Akk' matrix with two additional contributions 



- 2M sin(i^)„ A^,_fe To (14) 

+ t6 cos(i^)„A^,_fc+,A^, fg. (15) 

Eq. (14) corresponds to the usual term —J • A, where 

J„ =ak„Ck = 2t(5sink„ (16) 

is the quasiparticle current. Notice that in the presence 
of HF corrections to the quasiparticle dispersion (13) the 
quasiparticle current J^ is different from the quasiparticle 
velocity 
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In the usual Landau EL language, these two quanti- 
ties are related by the Landau FL i^f corrections as 



J a = Va{l + F^ /3) for an isotropic system in three 
dimensions'^. As we shall see, in our approach the anal- 
ogous role of Landau FL corrections will be played by 
the HS fields related to p-h fluctuations. Neglecting the 
diamagnetic contributions, that are not relevant for the 
following discussion, the leading terms in A of the ac- 
tion S'gl(A,T^-5) {T^^ = (j)" , (j)" , A) can be written in 
a compact form as: 

Sgl{A, T^^) = -ix^^ A^A^, + A^F"(T«S) (18) 

where Xjj is the mean-field current-current correlation 
function 



t-f(9) 



^^ ^ *^n + ^k+a - ilk- 



(19) 

and F°'{T^^) is a function of the T'^'^ fields which de- 
scribes the connection of the electromagnetic potential 
to the HS fields. The current-current response func- 
tion A"^ (q) can be computed from the partition function 
Z[A] = JpT^'5e-'5«i(A,T"^) of i-j^g jj^o^gi g^g 
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(20) 
where xli'i'l) simply follows from the quadratic term of 
Eq. (18), while Sxjjiq) is the contribution coming from 
the fluctuating modes coupled to A, and depends on the 



explicit form of the F{T^^) function: 
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Sxt'iq) = (i^"(T^*)F'^(T^*)) 



(21) 



Starting from the current-current response A"'^(g) (20) 
the paraconductivity simply follows from the dynamic 
limit (q = 0,a; — > 0) after analytical continuation of the 
Matsubara frequency ujm to the real frequency uj 



d<j^[lmSxfnio,q^O)/uj]^^o, 



(22) 



while the fluctuations contribution to the diamagnetism 
is connected instead to the static limit (q — > 0, w = 0), 



Sxd^~[Sx-,{<l,^ = 0)/(l^] 
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where SxL is the transverse part of the fluctuation cor- 
rection to the current-current response function. In the 
absence of fluctuations in the p-h channel, only the pair- 
ing field is coupled to A so that F(A) ^ pA^ and one 
recovers the standard AL correction (see Eq. (27) below) . 
In our case, one immediately sees from Eq. (9) and (14) 
that the 0^ field appears in the actions with the same 
structure of the electromagnetic potential A", i.e. it is 
coupled to the fermionic current. Thus, we expect that 
it will contribute to the Sxjjiq) correction (20) above. 

With lengthy but straightforward calculations one ob- 
tains that the effective action (18) at leading order in the 
gauge and HS fields is given by 
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where summation over repeated indexes is implicit and 
only the terms relevant for the following discussion are 
included. The c and c' terms in Eq. (24) follow from the 
expansion to leading order in q of the fermionic bubbles 
associated to the A^ term and to the AA^ and (p'^A'^ 
terms, respectively (see Fig. 1). In the former case for 
example one expands 



T 



^('i) = -MT.^di^^^^+q/2G- 



k+q/2 



H(w) - cq2 



(25) 
where the quasiparticle Green's function, Qf; = (iujn — 
^k)~^ J contains the full dispersion ^k- As a consequence, 
c is proportional to the second-order derivative of Q, 



which in turn scales as (d^k^) 



In the case of 



the c'q term instead one carries out a single derivative 
of the fermionic bubble which contains already a current 
insertion J, associated to each A or 0" field, see Eq.s 
(9) and (14) and Fig. 1. Thus, we have in a short-hand 
notation: 



c^id^kS^v 
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c'cxaa„9eL ^^'fJf (26) 



where vp, Jp are the average values on the Fermi sur- 
face. If we do not consider the interactions in the p-h 
channel leading to the HE correction of the band dis- 
persion we would get ^k = Ck (see Eq. (13)) and the 
velocity coincides with the current, so that c ~ c' . In 
this case, as mentioned above, the electromagnetic po- 
tential A is coupled only with the pairing field A via 
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where 



and 
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l/9-Xn{<l)/{'2tSy 



(29) 



c{q)=c'{l + Z{q)). 



(30) 



(b) 




FIG. 1: Fermionic bubbles: (a) is the Ti-{q) bubble associated 
to the A^ term, while (b) is the fermionic bubble relative 
both to AA^ and <j>"A'^ terms. The solid lines are the Green's 
functions Q containing the full quasiparticle dispersion ^. The 
dashed lines are the pairing field A. The wavy line represents 
either the electromagnetic potential A or the current field (f)" , 
both associated to the quasiparticle current J oc 9^". 



a term ^ cApA^, so that the fluctuation correction to 
the current-current response function (20) is given by the 
usual AL contribution^ (see Fig. 2): 

5x"^(g) =c2T^(2p + q)„(2p + q)^L(p)L(p + q) (27) 



where L{p) — (A* Ap) is the propagator of the SCF. Since 
c oc Wp is constant, the SCF contribution to the diamag- 
netism (23) and conductivity (22) is in both cases pro- 
portional to Vp , and the ratio Sa/6xd is expected to be of 
order 0(1). Such a result changes in the presence of p-h 



We notice that in Eq.s (28)-(30) we used a simplified 
notation valid in the g — s- limit relevant for the fol- 
lowing discussion, so that the current-current correlation 
function Xjj refers to its diagonal part x?f only. At fi- 
nite q the above expressions must be properly extended 
to account for the full matrix structure of the response 
functions. From the definition (28) and using Eq. (20) 
we compute the SCF contribution to the current-current 
response function. We find a correction to the mean-field 
value current-current response function (19) that leads to 
the RPA resummation of Xjj as 
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(<?) ^Xn{q) + Z{q)X]]{q) = 
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(31) 
and a second term that represents the SCF contribution, 
which generalizes the standard AL result (27) with a mo- 
mentum and frequency dependent vertex c{q): 

Sx"^iq) = T^c(9)2(2p + q)„(2p + q)f,L{p)L{p + q). 

p 

(32) 
From Eq. (19) one can easily see that in the dynamic 
limit Xjj = 0, so that Z = in Eq. (29) and from Eq. 
(30) c{q) — c' . In the opposite static limit instead we can 
rewrite Xjj in Eq. (19) as 



^^^ = -]^E(^k.^£)^a^/(/3a) 

k 

= |E-(Ce^)/(/3ek) 

k 



(33) 



FIG. 2: (color online) AL contribution of the SCF to the 
current-current correlation function. The dashed lines repre- 
sent the SCF propagator, while the red dots reduce to the 
constant c in the ordinary case Eq. (27), and to the momen- 
tum and frequency dependent vertex c[q) in the presence of 
current-current interactions, see Eq. (32). 

interactions. The elimination via gaussian integration of 
the current field (f)'^ from Eq. (24) leads to 



5gl(A,A) = --xj/(l + Z(g))A^A^,+ 



+ (g-i-n(u;)-fcq2)A;A 
- c(g)qA,,(A;A,_, 



By direct comparison with the self-consistent equation 
(12) for the density field one immediately sees that Xjj — 
2t5{JF/vF){(l)%lg)- From Eq. (29) it then follows that 
1 + Z = vf/Jf, leading to c{q) = c in Eq. (30). We then 
find that the vertices relative to the SCF contribution 
to the conductivity and diamagnetism are quantitatively 
different: 

c{q) ex vfJf '-^ i^tS + (j)^o)2St (q = 0,w^O) (34) 
c{q) ex vl r^ {2td + (j)^)^ (q^O,w = 0). (35) 

The difference in the two limits reflects in a difference 
in the overall prefactors of the SCF contribution to the 
conductivity and diamagnetism. Indeed, from Eqs. (22)- 
(23) one has that 
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i.e. the fluctuation conductivity is suppressed by the 
proximity to the Mott insulator by a factor that depends 
on the doping. 

We notice that in the presence of HF corrections the 
gauge-invariant form of the GL functional for SCF is re- 
covered in a non trivial way. Indeed, in the usual case 
the coupling to the gauge field A can be obtained by the 
minimal-coupling substitution q — > q — 2 A in the q^A^ 
term of the Gaussian propagator. This leads immediately 
to the term linear in A, cA • qA^, needed to compute 
the AL correction, see Eq.s (18)-(21). In our case by 
direct inspection of Eq. (24) one finds instead two differ- 
ent coefficients c,d in the q^A^ and A • qA^ terms, as 
we explained above. However, by integrating out the (^'^ 
field the coupling of SCF to the gauge field is described 
in general by a term c((7)A • qA^, (see Eq. (28)), where 
c{q) is given by Eq. (30). As a consequence, one recovers 
the minimal-coupling prescription only in the static limit 
(35) where c(g) — (?^ . 

According to the discussion below Eq. (1) above, the 
result (36) can be recast as an estimate of the ratio be- 
tween the correlation lengths extracted experimentally 
from paraconductivity and diamagnetism. For example 
using data reported in Ref. [13] one has that at T ^ 24 K, 
6a - 10^(17 m)-\ and bxd ~ 60 A/mT. Using in Eq. (1) 
^0 ~ 1 nm and d ^ \h A as appropriates for cuprates, we 
obtain that ksTld^l = 0.053 A/mT and e^ jX^MH =■ 
10'*(r2m)-i, so that Saj^xd oc Ca/^L ^ lO^^. Such a 
strong suppression of the paraconductivity with respect 
to diamagnetism, that has been rephrased in Ref. [13] in 
terms of the vortex diffusion constant valid only in the 
case of KT fluctuations, can be more generally attributed 
from Eq. (36) to the overall S^ factor due to the proxim- 
ity to the Mott-insulating phase. We checked that the 
estimate of the ratio (36) within the mean-field solution 
of the t — J model is quantitatively larger by a factor 
ten than the one experimentally found, since a prefac- 
tor ~ 2i/(/)g partly compensate the 5"^ suppression. Such 
a quantitative disagreement is reminiscent of analogous 
limitations of the mean-field approach already discussed 
in the literature in the contexts of other physical quanti- 
ties, as for example the scaling of the superfiuid-density 



depletion with doping^^. Moreover, at low doping^^'^^'^^ 
one should also include the boson fiuctuations neglected 
so far, which give a temperature Tb for the boson con- 
densation smaller than the one for the gap opening, lead- 
ing to a suppression of the critical temperature Tc with 
respect to its mean- field value Tmf- In such a regime, 
the static limit of the I + Z correction in Eq. (30) above 
could not be simply given by vp / Jp: nonetheless, the dif- 
ference between the static and dynamic limit still holds, 
possibly leading only to a quantitative difference with 
respect to the result (36). Finally, we notice that the 
experimental estimate of the SCF contribution to para- 
conductivity reported in Ref. ]13] is based on the scaling 
of the conductivity at frequencies large enough with re- 
spect to quasiparticle dissipation. As a consequence, it 
is worth making a comparison with the present results, 
that have been derived in the clean limit. 

In summary, we analyzed the SCF contribution to con- 
ductivity and diamagnetism in the presence of current- 
current interactions, by using as a paradigmatic exam- 
ple the slave-boson formulation for the t — J model. 
By explicitly constructing the GL fluctuation func- 
tional in the presence of HF corrections we showed that 
current-current interactions, needed to recover the gauge- 
invariant form of the GL functional, modify the transport 
coefficients leading to a momentum and frequency depen- 
dence of the vertex c{q) entering the AL expression for 
the SCF contribution. Since different limits are involved 
in the definition of paraconductivity and diamagnetism, 
we obtain a different prefactor in the two cases, with a 
suppression of the paraconductivity because of the prox- 
imity to the Mott-insulating phase, as recently shown 
experimentally in cuprates^'^. Even though a mean-field 
approach to the t — J model is not satisfactory from the 
quantitative point of view, the different strength of SCF 
contribution to conductivity and diamagnetism is more 
general, since it is a consequence of the existence of a size- 
ble difference between quasiparticle current and velocity. 
A quantitative comparison with experiments remains an 
interesting theoretical challenge, that certainly deserves 
further investigation. 
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